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We study the charge character of the Andreev bound states ( ABSs) in one-dimensional topological 
superconductors with spatial inversion symmetry (SIS) breaking. Despite the absence of the SIS, we 
show a hidden symmetry for the Bogoliubov de Gennes equations around Fermi points in addition 
to the particle-hole symmetry. This hidden symmetry protects that the charge of the ABSs is solely 
dependent on the corresponding Fermi velocities. On the other hand, if the SIS is present, the ABSs 
are charge neutral, similar to Majorana fermions. We demonstrate that the charge of the ABSs can 
be experimentally measured in the tunneling transport spectroscopy from the resonant differential 
tunneling conductance. 
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Introduction - In superconductors (SCs), Bogoliubov 
quasiparticle (BQ) is a coherent superposition of the elec- 
tron and hole; hence the charge carried by a BQ is not 
conserved [1]. The effective charge of the BQ can be 
identified as the expectation value of the charge operator. 
Based on this definition, the charge of the quasiparticle 
excitation above SC gap is generically a continuous func- 
tion of the energy and SC order parameter, and therefore 
it is expected to be sensitive to perturbations [1]. 

A peculiar subgap quasiparticle state in topological 
SCs is Majorana fermion (MF), which is charge neu- 
tral and equal to its own antiparicle. Driven by the 
exotic properties that MFs obey non-Abelian statistics 
and have potential applications in fault tolerant topolog- 
ical quantum computation i2]-t4j, the search for MFs in 
realistic condensed matter systems has become an excit- 
ing pursuit. Candidates shown to exhibit topological su- 
perconductivity and host MFs include two dimensional 
(2D) p + ip IMS] and ID p^-wave SCs [III [H], topo- 
logical insulator/SC heterostructures [T3'-fro, and spin- 
orbit (SO) coupled semiconductor/SC heterostructures 
[T71 - f25| . MFs are self-Hermitian and therefore charge 
neutral, which in SC is guaranteed by the particle-hole 
symmetry [5]. Since the self-Hermitian property exists 
only for zero modes, a natural question arises: what is the 
charge of a generic (non-Majorana) Andreev bound state 
(ABS)? Furthermore, in the tunneling transport spec- 
troscopy it is predicted that a single MF can induce the 
differential tunneling conductance (DTC) peak of height 
2e^/h at the resonant Andreev reflection which 
is a consequence of the self-Hermitian and charge neu- 
trality of MFs. Is such peak value of DTC unique to the 
MFs? These fundamental questions motivate us to look 
more closely at the ABSs in SCs. 

In this Letter we focus on the ABSs at Josephson junc- 
tion of ID topological SCs. We show a hidden symmetry 
for the Bogoliubov-De Gennes (BdG) equations around 
Fermi energy, by which we predict that the charge of 
ABSs solely depends on the Fermi velocities, regardless 



of other details. We also demonstrate that the charge can 
be measured by DTC at the resonant Andreev reflection. 

Generic theory - We start from the generic mean- 
field Hamiltonian for ID spinless p-wave SC O fT2] : 
^ — — M)cJ,Cfc-|- J dxA{x)c{x)c{x + a) + h.c., where 

A is the induced SC order parameter, fi is the chemical 
potential, a is the lattice constant, c and are annihi- 
lation and creation operators of electrons, respectively. 
A broken spatial inversion symmetry (SIS) is generically 
allowed for the present study and is characterized by the 
single electron dispersion relation £k ^ £-k, which leads 
to a discrepancy between Fermi velocities VfL and Vffj 
for the left and right movers. Different mechanisms for 
SIS breaking will be studied later in lattice model and in 
semiconductor nanowire systems. 

In the parameter regime that |A| <C Ep with Ep the 
Fermi energy measured from the bottom of the band, 
we can expand the Hamiltonian around the Fermi points 
by the transformation c(x) = Cij(a:)e*'^^^ -I- CL(a;)e~*'^^^, 
and keep the terms up to the first order of l/{kii,LO 
with ^ the coherence length of the SC. Here ki, and 
fc/j are respectively Fermi momenta at left and right 
Fermi points. With the help of the Nambu bases ^+ = 
[c-r{x),c^^{x)Y' and -0- — ['^l(x) , c^-pi{x)Y' , we rewrite 
the Hamiltonian in explicit particle-hole symmetric form 
= / dxi;X{x)U+{x)i!+{x) + J dx^plix)'H-ix)^p_{x), 
with ^{-{x) related to H+{x) by particle-hole transfor- 
mation 'H-(x) = —Ty'H*^_{x)Ty, and 



H+{x) 



-ivfRdx iA{x) 
-iA*{x) ivfLdx 



(1) 



The explicit form of A(x) depends on the electron dis- 
persion relation S^. and Fermi momenta, but as shown 
below, it does not affect the charge of the ABSs. 

The ABSs can be studied by considering a Joseph- 
son junction formed around a; = 0. We denote by 
the ABS wave functions ^,y±{x) = [uiy±{x),Vi^±{x)]'^ 
for the BdG Hamiltonians H± with eigenvalues £^± 
[v — ... — 1,0,1...). The Bogoliubov operators of the 
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ABSs are then defined as b^^ 
v^±(x)c'j^^j^{x)] and the electric charges are calculated by 
e*_|_ = e J c?x[|M^-|-(a;)p — |wi/±(a;)p]. However, the wave 
functions of ABSs are not analytically solvable in the 
generic situation. Fortunately, without solving the wave 
functions, we shall show that u^± (x) and v^± (x) sat- 
isfy a universal relation for all ABSs: \u^^{x)/vi,-^-{x)\ — 
\u^^{x)/v^_{x)\^^ = (w/i/w/fl)^/^. We then have 



s*± = el = ±e{— —). 



(2) 



We prove the above result with the BdG Hamiltonian 
Hj^. The SC order parameter is written as A(a::) = 
Ai(a;) + iA2{x). Our main goal is to show that the re- 
lation VfR\ui,+ {x)\'^ = VfL\vu+{x)\'^ is valid for all ABSs 
in generic case. For this purpose we define Gi,+ {x) — 
\uij+{x)/viy^{x)\'^ , and derive from BdG equations that 



dGt,+ {x) Av,,+ {x)\'^ Im(Au*+t;,,+ ) 



dx 



VfR 



VfL 



(3) 



Note that Im(2Au;^_^ 

bound states. Then if one has 
only solution is \u^+{x)/v^+{x)\'^ 



) = VfJid\u^-^-\'^/dx 7^ for 

= 0, the 
VfL/vfR. On the 
0, where 



iu+{X)/ ui,+ y 

other hand, we notice that 92;/,y+(Ai, A2) 
/i.+ (Ai,A2) = VfRlu^+l"^ - VfL\v^+\'^. Therefore the 
function f^^{Ai, A2) is independent of position. To de- 
termine the value of /,y+(Ai, A2) we examine the asymp- 
totic behavior of the ABS wave function. We take the 
boundary condition that the order parameter A(x — > 
±(X)) is finite. The asymptotic behavior of u^^{x) and 
v^+{x) can be generically described by M^+(a; +00) ^ 
^g-7i^-H*7i^ and v^+{x +00) ~ ^e-''2^+*^2x^ j.g_ 

spectively. By examining the BdG equations of u,y-^(x) 
and v^+{x) at a; -l-oo one can verify that the expo- 
nents must satisfy the relations 61 = 62 and 71 =72- 
This confirms that G^-^{x — > +00) is a constant; hence 
.+ (a^^+oo) ^ Q According to Eq. ^ it follows then 
Vfh/vfR. Thus we get /^+(Ai, A2) = 0, 



dx 

G^+(x — ^ 00) 

which is valid for the entire position space. Finally we 



reach the result VfR\ui,+ {x)Y 



VfL\Vv+ 



(x)| , complet- 



ing the proof. Eq. ^ gives the charge carried by any 
ABS in the generic case. We emphasize that this proof is 
valid for any bound state in the bulk region (with wave 
function exponentially vanishing for x — ^ ±00). 

The results of the charges e*^ indicate a hidden 
symmetry for the BdG Hamiltonians. Let — 
[uv±{x),v^±{x)Y satisfy 7^±$y± = £^±<^^±. We denote 



by /C and Us^^^ = 



)X 



the complex conju- 



gate and the U(l) unitary transformations, respectively. 
Furthermore, we define the Lorentz boosts C± — e^^'^' 
to rescale u^±{x) and Vy±{x), with A± = ± ln(wL/wi?,)"'^/^ 
and T, 
It 



x^y^z the Pauli matrices acting on Nambu space, 
can be verified that the Hamiltonians are invari- 



ant under transformations Hi/±'H±Hj^_^ = "H^, where 



via $i/± 



U£^,JCtxJC±, while the wave functions transform 



One can show that the bound 



states of y.^ {T~LJ) are non-degenerate. This leads to 
— e*''±4>*_|_, with •q± arbitrary constants, which im- 
plies that 



Ui,±{x) — e^^v 



)x 



(4) 



This formula recovers the charges e^. in Eq. ([2| . Here we 
have shown a remarkable result that the electric charge 
of an ABS is solely determined by the SIS breaking term 
which is quantified by the discrepancy of the Fermi ve- 
locities Vf^ but independent of all other details. This 
result is protected by the hidden symmetry of the BdG 
Hamiltonians 'H± under the transformations Hj,±, and is 
valid when l/(fcL/i^^) ^ 1. Particularly, for the inver- 
sion symmetric ID SC with VfR = VfL, the ABSs are 
charge neutral, regardless of the details of the SC. 

Lattice Model for spinless p-wave SC - Next we study 
concrete systems with SIS breaking. We consider first 
the lattice model for ID topological SC described by 



H ^ - Y^ tic\cj - t2e''^^'clcj -fil^c^c. 



(5) 



where ii is the nearest hopping, and the complex next- 
nearest-neighbor (NNN) hopping i2e'^'^ (with 9i^i±2 — 
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FIG. 1: (Color online) Numerical (solid lines) and analytic 
(dash-dotted lines) results for the charge (in unit of e) of the 
ABSs versus 6^0 (a), SC gap |A| (b), and Josephson junction 
phase difference (l>o (c), with ti rescaled to be dimensionless 
and taken as ti — 10. (d) Energy of the ABSs versus (f)o, with 
6*0 ~ 2n/5 and |A| = 1.0. Other parameters are taken as /i < 
0, t2 = 0,0.5, 1.0, 2.0, ?io = tt/2 for (a), t2 = 0, 1.0, 2.0, (?io = 
7r/2,6lo = 37r/2 for (b), and = 0,2.0,6lo = 27r/5, |A| = 1.0 
for (c-d). 
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±00 ) breaks SIS, as in the Haldane model [57| and Kane- 
Mele model [28^. Recent proposals show that the com- 
plex NNN hopping with appreciable magnitude may be 
obtained in the silicene [IHl |30] , and realized by doping 
localized spins in graphene [3T] . The Hamiltonian Eq. ^ 
becomes Kitaev model when ^2 = [T^], which can be 
effectively achieved by e.g., s-wave SC proximity effect 
in the ID semiconductor nanowire (NW) with Rashba 
SO interaction [T5H11]- Transforming H into k space 
yields H = ^^{£k - + I]fe[^sin(A:a)cfeC_fc + h.c] 

with Sk = —2ti cos(fca) — 2t2 cos(2/ca + 6*0). The electron 
dispersion relation £k ^ £-k (for /c ^ 0) possesses no 
inversion symmetry when do 7^ n-n and |i2| > 0. 

We present in Fig. [l] the analytic and numerical re- 
sults for the ABSs localized on a Josephson junction with 
phase difference (pQ. Fig.[T](a-b) shows that the numerical 
solution is well consistent with Eq. ([2| in the parameter 
regime that \/^\/Ef < 1/10 (equivalent to the condition 
^/{kh^RO < 1/10). Fig.jl (c-d) confirms that the charge 
is independent of ^0 and the energies of ABSs. It is 
noteworthy that the crossing point of the ABS spectrum 
{v — 0) is shifted away from (jjQ = n when SIS is broken 
(blue stars in Fig. [I] (d)), while the zero energy crossing 
itself is topologically stable, at which point two zero MF 
modes can be introduced [HI [19] . 

Semiconductor Nanowire - For the semiconductor 
NW/SC heterostructure, the SIS is broken when the Zee- 
man field {Vx,Vy) = Vb(cos6'o,sin6'o) has a tilt angle (6'o) 
relative to the NW direction. The Hamiltonian reads 



V„/E =3.0,4=11/2 



|A 1/E =0.3, V„/E =3.0 



H 



E 



dxcl,{3 



-3. 

-2m* 



H + {iXB.dx + Vy)(Ty (6) 



{x) -\- I dx[/S.s{x)c^[x)c\^{x) + h.c. 



where m*,Xji and are effective mass of electrons, 
spin-orbit coupling coefficient, and induced s-wave SC or- 
der parameter, respectively. Note that under the spatial 
inversion transformation x —x, one has {kx,cry) — >■ 
{—kx,—cfy), and Gx cFx- Therefore the term VyUy 
breaks SIS and results in a discrepancy between Fermi 
velocities VfL and w/_r. For the case 9q = 0, the above 
system is isomorphic to Kitaev model and support Ma- 
jorana end modes when Vx > |Asp -I- /x^ [TU1I5T] . 

For a large Zeeman field relative to |As| and one 
can project the Hamiltonian to the lower subband and 
linearize the electron dispersion relation around Fermi 
points. The charges are then obtained according to 
Eq. ([2|. The results are confirmed in the numerical cal- 
culation given in Fig. |2] (a-b) . 

Tunneling transport spectroscopy - Now we propose 
to detect the ABSs with the tunneling transport spec- 
troscopy. A single normal metallic (NM) lead with volt- 
age eV is coupled to the ABSs at the Josephson junction 
(Fig. [3]) in the tunneling regime. The coupling Hamil- 
tonian is described by Ht = / dxt* {x)d^ {x)c{x) -t- /i.e.. 





FIG. 2: (Color online) Numerical (solid curves) and analytic 
(dashed lines) results for the charge (in unit of e) versus tilt 
angle 60 (for (a)) and Josephson junction phase difference 4>(i 
(for (b)). The spin-orbit coupling energy Eso = rn*\\ = 1, 
and the chemical potential = 0. 



where d{x) is the electron annihilation operator in NM 
lead, and the tunneling coefficient \t{x)\ > in the region 
|a;| < dn/2 with dn the width of the NM lead/SC contact 
region. Projecting the operator c{x) onto the manifold 
of ABSs yields the effective tunneling Hamiltonian 



(7) 



where f^^^ = / dxtl{x)ul±{x)e ''^"/^^ and gk,„± = 
/ dxtk{x)vl^{x)e''''^/'^^, with tk{x) = t{x)e^''^ . The 
tunneling current is calculated by the rate of change 
of the electron number N = d\dk in the NM 
lead: / = —eN — —^[Ht,N], which can be stud- 
ied with the Keldysh formalism. Then we define 
the Keldysh contour Green's functions for ABSs by 
Qmp'(^,^') = -^{TK[b^,iT)bl,{T')]), with the index ^l = 
(;/, ±), and the free lead electron Green's functions 
by GO(r,r') = -z(TK[4(r)4.(T')])o and GO(r,T') - 
—i{TK[dj.{T)dk{T')])Q. Here {■)o represents the situation 
with Ht = 0. The current is obtained by 

-|-|/^Tr[(E('')Q-QS('^)):], (8) 

where the trace is performed in the space spanned by 
modes. The self-energy matrices in real time space 



M 




4>o 

(b) 



FIG. 3: (Color online) Sketch of the tunneling charge trans- 
port between the NM lead and the Josephson junction. A 
wide contact (a) and a narrow point contact (b) between the 
NM lead and the SC are considered. 
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Ekh,,ft,'Gl{r,r') and ^l^^,Ar,r') 



Efc 9k fj.dk.fj.'G^i''': '^')- We derive the equation of motion 
for Qt,^'{T, t') and obtain {idr - E^^'^s - S)Q = 1, where 
S = S^''' + S^''^ and E'l'^s = diag{ . . . , , ...}. The solu- 
tion reads Q = + Q"SQ, with Q" = (w - E'^''^s)-i. 

For the NM lead, we consider the wide band limit 
that the transition matrix elements /fc,^' and gk.^i' are 
weakly energy dependent [32]. In this case the self- 
energies are purely imaginary and we denote the retarded 
components by S^''jfi(a;) = |T^^,,i(w) and S^'^l^2(^) = 



2 A'A''.2 



determine 



MM 

(-w). Physically T^^.^ and T^^/^2 
the tunneling transmission probabilities for electrons and 
holes, respectively. 

It is interesting to consider two different geometries of 
the NM lead characterized by the width (i„. First, we 
study the wide contact regime (Fig. [3] (a)) with the width 
dn greater than the localization length ^abs of ABSs. 
Note that the Fermi wavelength in NM lead is much less 
than Iabs- This leads to a fast oscillating behavior for 
the phases of /fe^^ (5fc,/i) versus k [25]. The off-diagonal 
elements of the self energy then vanishes, and T^^/ i w 
'^MM' Sfc /'^2;da;'t^(a;)ife(x')u* (x)u^(x')J(efc - w) and 
T^^',2 ~ (^MM'Sfe/ dxdx'tl.{x)tk{x')vl_{x')v^{x)5(ek-i^): 
with efc the free electron energy in the NM lead. 
The retarded Green's function for ABSs is given by 
{QVM = + with = (T^^,i + T^^,2)/2. 

Bearing these results in mind, we obtain the DTC at 
zero temperature 



dl 
dV 



2e2 



E 



.iT 



{eV- 



y2 ' 



(9) 



One can verify with Eq. Q that T^^^i/T^^^2 — 7^- Let 
i£'min be the minimum energy spacing for the ABSs, which 
depends on 0o (Fig. [Jd)). For < S^.^, the DTC 
has peaks at eVm = Ct^, corresponding to the resonant 
Andreev reflection induced by the ABSs. The resonant 
DTC for the peak values takes a simple form 



dl_ 
dV 



2e2 
h 



2ef 



(10) 



which measures the charge e^. carried by ABSs. 

Second, for the narrow point contact regime 
with dn ^ 'abS: the coupling can be approx- 
imated as t{x) ~ tQdnS{x), which follows then 
T^A^M « u;iO)u^,iO)\to\^dlN{tu) and T^^,,2 » 
v*^{0)Vf,'{0)\to\'^dlN{u}), with N{uj) the density of states 
in NM lead. The main difference in the present regime is 
that the off-diagonal elements T^^/ (^ ^ are gener- 
ically nonzero, which may lead to interference effect in 
the tunneling transports through different ABSs. Nev- 
ertheless, when T^^, <^ ^mim the off-diagonal terms are 
not important for the resonant Andreev reflection and we 
show that peak values of the DTC have the same form 



generic formula of the DTC through resonant Andreev 
reflection induced by each ABS with effective charge e*. 



dl/dV(2e^/h) 



dl/dV(2e7h) 



(a) 




(b) 




FIG. 4: (Color online) DTC dl/dV (in units of 2e^/h) versus 
the bias eV (in units of To). The blue and red curves respec- 
tively correspond to the inversion symmetric case {vfL ~ vfR) 
and the case of Vfi^/vfR = 0.44 {t2 = 0.2ii and So = 7r/2 in 
lattice model), (a) The Josephson junction of a short length 
il < i)/l/|A|) is considered, with only two ABSs localized on 
the junction; (b) There are six ABSs locaHzed on a long junc- 
tion (Z ~ 3«/r/|A|). The equal height of peaks tells that the 
absolute value of the charge are the same for different ABSs 
in each parameter regime. The black dotted lines correspond 
to {dI/dV)m = 2(e2 - e*^)/h. 



Fig. [4] (a-b) shows the DTC as a function of bias eV. 
The peak values of the DTC at eVm = £fj, confirm the 



results in Eq. 10 for both the wide contact regime (solid 



as in wide contact regime. Therefore Eq. ( 10 1 gives the 



curves) and narrow point contact regime (dashed curves) . 
In particular, for Vfu = VfL (blue curves) one has Cj. = 
and the peak value of the DTC at resonant Andreev re- 
flection is given by 2e^//i. Physically, this can be under- 
stood that an incident electron leads to the unity reflec- 
tion of a hole at the resonant Andreev reflection induced 
by a neutral ABS. This process is the same as the MF 
induced resonant Andreev reflection (MIRAR) [21] . Our 
results show that the induced peak value 2e^/h of the 
DTC is universal for any single neutral ABS in ID SCs. 

In summary, we have presented a profound relation be- 
tween the charge of ABSs and SIS in ID topological SCs 
with small SC gap. A hidden symmetry is revealed for 
BdG equations, which protects that the charge of ABSs 
is solely determined by the Fermi velocities, regardless 
of other details. The charge of the ABSs can be mea- 
sured by tunneling transport spectroscopy. Particularly, 
for system with SIS, the ABSs are charge neutral and 
each induces a DTC peak of height 2e^ /h at the reso- 
nant Andreev reflection with finite bias, similar to the 
cases for multiple MFs [24ll25]. 

We thank K. Sun, R. M. Lutchyn, T. D. Stanescu, K. 
T. Law, V. Yakovenko, Z. X. Liu, and X. Liu for helpful 
discussions. This work is supported by funding from JQI- 
NSF-PFC, Microsoft-Q, DARPA, and QUEST. 
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Andreev Bound states in One Dimensional Superconductor — Supplementary Material 

In this supplementary material we provide the details of some results in the main text. 

Proof of the charge carried by Andreev bound states 

The BdG equations of u^+{x) and v^+{x) are given by -^u^+{x) — — i^j^t;,y+(x) + and ^Wy+(a;) — 

i^^^j^Uiy^{x) — i^^v^+{x), from which we get further their complex conjugate counterparts by -^u*_^^{x) — 
— i^^ul^{x) and ji^vt+i^) = + ^ITJ^'^t+i^)- With these formulas we can derive that 

dG^+{x) ^ Av^+{x)\^ Im(AM*^t;^+) ^^^^ 

dx VfR VfL \v\'^ 

It can also be verified from the BdG equations that ^(w/_r|u^_|_P — VfL\v^+\'^) = 0, which suggests that 

/,+ (Ai, A2) = P - w/ik.+ P (12) 

is a function independent of the position. The value of /(Ai,A2) can be determined by examining the asymptotic 
behavior of the ABS wave function. Note that the asymptotic behavior of u^^ix) and v^^ix) can be described by 
— > +00) ^ g^j^j Vyj^{x — > +00) ^ A^^-ii^+^iix ^ respectively. Here ^1,2,^1,2,71,2, and 71,2 are 

constants. For x — >■ +cxi we take that A(a; — )• +00) — Aq is finite, and from BdG equations we have 

where the higher-order infinitesimal terms are neglected, j 7^ j' = 1,2, mi — w/i?,, and zu2 — VfL- Comparing 
both sides of the above equations we observe that the exponents must satisfy the relations Si ~ S 2 and 71 =72- 
This confirms that G^+{x — +00) is a constant, and thus '^^•'+'^^^+°°'> = q According to Eq. |ll| it follows that 



Gy+{x — )■ 00) — \u^+{x)/v,j+{x)\ — VfL/vfR. Together with Eq. (12| we have /^4.(Ai, A2) = 0, which is valid for 



the entire position space and finally we reach = vji^\vi,^{x)\'^, completing the proof. 

The hidden symmetry of the BdG Hamiltonians 'H± 
A. Non-degeneracy of the bound states for %+ 

Before turning to the hidden symmetry, we show a basic property that the bound states of T-L+ (the same for ?^_) 
are non-degenerate. Let $1 = [ui{x)^vi{x)Y' and $2 = [^^2(2^), ^2(0;)]^ be two degenerate bound states of T-L+ with 
energy 8. Then we have j^Uj{x) = -i^Vj{x) + i:^Uj{x) and £vj{x) = i^Uj{x) - i:;^Vjix) with j = 1,2. 
From the BdG equations one can show that 

-^{uiV2 - U2V1) = —){uiV2 - U2V1), (14) 

dx VfR VfL 

which has the solution uif2 — U2'yi = Ce "f^ "f^ with C a constant. Since $1.2(3;) are bound states, we have 
M1W2 — U2V1 — > for a; — > 00, which implies that C = 0. Then one gets 

Vi{x) V2{x)' 

Namely, $i(x) and $2(2;) are the same state, and therefore the bound state spectrum of is non-degenerate. 
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B. Hidden symmetry 

We consider 71+ first. Let $y+ — [u„+{x) , v^-^-{x)]'^ satisfy H+'J>,y+ = £u+'^v+- First, we apply tlie transformation 



where lAs^ 



a>,+ = (W£„^/c)$:+, (16) 

,^ — ^ '/« "/J^ ''^ is the U(l) transformation and /C is the complex conjugate operator. This leads to 

(17) 



VfR VfL 



dx ^ VfL ^ VfR 



(18) 



Furthermore, we define the operator = Tx^+, where £+ = e'^+'^* is a Lorentz boost to rescale u^+{x) and v,j+{x), 
with A+ = hi{vfL/vfiiy^^, and t^; ^ ^ are the Pauli matrices acting on Nambu space. The operator V-^. takes the form 





e-^+ 



Under the transformation 



which gives = [■D*^(x), {t*_|_(x)]-^, we obtain 



d 

dx 



w*+(x) 



.A (a;) , s .£u+ I* / s 
-I u^^{x) + I ^^+(2;), 

VfR ^ VfR 



(19) 



(20) 



(21) 



— u^+(a:) = I v^+{x)~i u^+{x). 

dx VfL VfL 



(22) 



The above equations can be rewritten as = E^+^l^. Therefore is still the eigenstate of 7i+ with the 

energy Note the bound states of are non-degenerate. This leads to = e"'+$*_|_, with 77+ arbitrary 

constant, which gives 

u,+ {x) = e^+w:+(a;)e'''e'^"+^^"^'". (23) 

Similarly, for the Hamiltonian H_(a;), we introduce the similar transformations V-jK, and Us^_, where Us^_ — 
Q^''~^^SR "fL -'^ and V- = t^C- with £_ = e^-'^' and A_ = — ln(w/i/w/fl)^/^. Under these transformations we also 
find 'H-^l_ — which leads to — e"'"$*_, and therefore 

u,^{x) = e^-t>:_(x)e"'-e'^''"'^"^^". (24) 
The above consecutive transformations for "Hi (x) can be summarized that under the transformation 

*.±(a:) = (Z^f„±/CT,£±)l:±(a;), (25) 

the Hamiltonians are invariant 



(26) 



Therefore from Eqs. (23 24| we have j^Jf^^fyjr = e^^+ = ^ and \llZ{x)\^ = e^^" = ^ff ■ This explains why \u^±{x)\^ 
and |wi/±(a;)P have such a simple relation and the charge only depends on the Fermi velocities for the linearized BdG 
Hamiltonians. 
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Differential tunneling conductance 

From the formula / = —^[Ht, N], we obtain the tumiehng cmTent by 



^ = ^ E E ifl,G<,iO, 0) - h,,G<^{0, 0) - gk^,g<,iO, 0) + glJ^^iO, 0)], (27) 

k fJ.=± 

where the mixed Green's functions are defined by Gkf_i{T, r') = — «(TR-[dfc(T)&j^(T')]), G^fc(r, r') — --i(TR-[6^(T)dj.(T')]), 

Gtj.k{T,T') — — i(T/f [&^(T)dfe(T')]), and Gkiiij.T') = — i(Ti<'[dfc(''")6j((T')]). In the first order approximation we obtain 
for the lesser Green's functions that 



and 



Gtkiry) ^Y.j dr"[h,,>Q,Ary)Gl{T",r')]'' 
Qt^r^) = - JlJ dr"[GUr,r")gl,,Q,',ir"y)y 



(28) 
(29) 

(30) 
(31) 



The tunneling current then recasts into 

^ = ^E/ dT[^i%,AO^r)Q,,^iT,0)-Q^^,{0,T)j:l%,{T,0) + 

+Q^,, (0, t)J:1%^^{t, 0) - S^'l^^CO, r)Q^,^{r, 0)] <, (32) 

where Y.^^^, j^{t,t') = Efc /fc.p/fe.p'G'°(T, t') and s|^^, 2(r,r') = 'Ek 9l,,gk,,,'Gl{T,T') are the corresponding self- 
energies. The Eq. (32 1 can also be written as I = Ii + I2, where 



e 



hi ^Tr[(Q«-Q:J)sS<+Q<(sS^-SS^)], (33) 



I2 = 



/gTr[(E^'')Q-QEri1 



hj 277 

I I £Tr[(Q^-Q^)s(t^<+Q<(s£)-^-S^t^«)]. (34) 



The Dyson equation of Qfj.f_i'{T,T') can be derived through idrQfifi'(T,T') ~ (^^^'(^(t — r') + i{TK{[H,b^]{T)b^^,(T'))), 



which follows that (idr - E^'^s - S)Q = 1. Here S = ' + S^^ and E^^'^e = diag{..., f^, ...} is the diagonal matrix 
composed of eigenvalues of the ABSs. The solution reads Q = Q° + Q°SQ, with Q° = {uj - E^^'^s)"!. 

For the NM lead, we consider the wide band limit that the transition matrix elements /fe,^' and gk.fj.' are 
weakly energy dependent [I]. In this case the self-energies are purely imaginary and the retarded components 
read J:^^I^}'^'{uj) = fT^^,_i(u;) and S^'),l^2(w) = IT^m'.sI-'^). where T^^/,i(w) = 27r /fc,^/*_^,5(w - efe) and 
'^nn',2i^) = '^'^J2k9k,fj.9kti'^i^ ~ ^k), with Efe the free electron energy in the NM lead. The lesser components 
= ^r^.^'A^)fi^ ~ eV), yI%{lo) = zT^^,,2(-w)[l - /(c. - eV)], and Q< - Q5S<Q^. This leads to 

/i = ^ 1 ^Tr{Q^(^)Y2(c.)Q^(^)Yi(c.)}[l-/(c.-el^)], (35) 

/2 = ^ 1 ^Tr{Q«(a;)Yi(w)Q^(c.)Y2(w)}[l-/(c.-eF)]. (36) 
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It can be verify that h = The differential conductance is then given by 

^ - ¥/ gT.iQ«,e^)r.Q^e.,r.,S<^. (37) 

The Fermi wavelength in the NM lead is much less than the ABS localization length Iabs, and also typically 
much less than k]^\ in nanowire systems. For the wide contact regime with the width d„ of NM lead greater 
than ZabS) the functions fk^^ {9k, ij.) exhibit a fast phase variation versus k, and the off-diagonal elements of the 
self energy vanishes. We then reach that T^^/.i « '^/j/j' X^fc / dxdx't*j^{x)tk{x')u*^{x)u^{x')5{€k — and T^^/_2 ~ 
(5^^' X]fc / dxdx' t*f.{x)tk{x')v*^{x')v iji,{x)5{ek — oj). The retarded Green's function for ABSs {Q^)l^^{oj) = w — + iT^, 
with = (T^^^i + T^^_2)/2. With these results we obtain the DTC at zero temperature 

rlT 9p2 

^ = ^TV{Q«(ey)T2Q^(ey)Ti} 

which is the Eq. (9) in the main text. For <C f^j^ with fmin the minimum energy spacing for the ABSs, the DTC 

has peaks at eVm ~ ±^/i, with the peak values given by (j^) = — -7^, which measures the charges e^. carried 
by ABSs. 
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